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Abstract 

We study the thermodynamics and phase structure of asymptotically flat non- 
dilatonic as well as dilatonic black branes in a cavity in arbitrary dimensions (D). 
We consider the canonical ensemble and so the charge inside the cavity and the 
temperature at the wall are fixed. We analyze the stability of the black brane 
equilibrium states and derive the phase structures. For the zero charge case we 
find an analog of Hawking-Page phase transition for these black branes in arbitrary 
dimensions. When the charge is non-zero, we find that below a critical value of the 
charge, the phase diagram has a line of first-order phase transition in a certain range 
of temperatures which ends up at a second order phase transition point (critical 
point) as the charge attains the critical value. We calculate the critical exponents 
at that critical point. Although our discussion is mainly concerned with the non- 
dilatonic branes, we show how it easily carries over to the dilatonic branes as well. 
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1 Introduction 



Over the years much is known about the thermodynamics and phase structure of black 
holes in asymptotically AdS space. The reason is that the black holes in AdS space, unlike 
those in flat space, are thermodynamically stable [1]. Further, more interests were drawn 
with the advent of AdS/CFT correspondence [2, 3, 4] as the AdS black holes provide a 
good laboratory for testing the correspondence at finite temperature [5]. So, for example, 
the AdS black holes are well-known to undergo a Hawking-Page phase transition [1] and 
by AdS / CFT this corresponds to the confinement-deconfinement phase transition in large 
N gauge theory [5] . Similarly, the phase structure of the charged AdS black holes which 
includes, in the canonical ensemble, a similarity with the van der Waals- Maxwell liquid-gas 
system, can also be understood from the dual field theories [6, 7]. 

However, it was pointed out [8, 9] that the above mentioned phase structure is not 
unique for the asymptotically AdS black holes. In fact, very similar phase structure was 
shown to arise for the suitably stabilized asymptotically flat as well as asymptotically dS 
black holes [8]. As, the higher dimensional theories like string or M-theory admits higher 
dimensional black objects like black p-branes [10, 11, 12], it is natural to ask what kind 
of phase structures do they give rise to - do they have a similar phase structure as the 
black holes or they have a different phase structure altogether? 

Motivated by this we looked into the thermodynamics and the phase structure of black 
p-brane solution of D-dimensional gravity coupled to (p+ l)-form gauge field. In the be- 
ginning, we consider only the non-dilatonic branes, (when D = 11, they correspond to 
M2 and M5 branes of M-theory and when D = 10, it is the D3-brane of string theory) 
and then towards the end we show how the analysis carries over to the dilatonic branes 
as well. The solutions we consider are as usual asymptotically flat and so they are ther- 
modynamically unstable and an isolated black brane would radiate energy in the form 
of Hawking radiation. In order to restore thermodynamic stability so that equilibrium 
thermodynamics and the phase structure can be studied, we must consider ensembles 
that include not only the branes under consideration but also their environment. As 
self-gravitating systems are spatially inhomogeneous, any specification of such ensembles 
requires not just thermodynamic quantities of interest but also the place at which they 
take the specified values. In other words, we place the brane in a cavity a la York [13] 
(see also [14, 15, 16, 17, 18]) and its extension in the charged case [19]. Concretely, we 
will keep the temperature fixed at the surface of the cavity and as the black p-branes are 
charged under the (p+ l)-form gauge field, we will keep the charge inside the cavity also 
fixed. This will define a canonical ensemble and we will study the phase structure of the 
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black p-branes in this ensemble. 

After some generalities we first consider the case when the charge in the cavity is 
fixed to be zero. In this case we find that there is a minimum temperature below which 
no black brane state exists in equilibrium inside the cavity. But above this temperature 
there exist two black brane states with different radii. The larger one is locally stable 
and the smaller one is unstable. The locally stable or 'supercooled' large black brane will 
eventually decay to energetically more favorable state 'hot flat space'. There is a phase 
transition temperature at which the 'hot flat space' and the corresponding large black 
brane can coexist 4 . But when the temperature rises above this transition temperature 
larger black brane becomes globally stable and the 'hot flat space' can decay to the black 
brane. The small black brane still remains unstable. As the temperature rises more, the 
size of the small black brane decreases and that of the large black brane increases and at 
infinite temperature the size of the small black brane goes to zero 5 whereas, the size of 
the large black brane coincides with the size of the cavity. This situation is analogous to 
the Hawking-Page transition of AdS black holes or asymptotically flat or dS black holes 
in a cavity [1, 8]. 

When the charge is non-zero but fixed inside the cavity, the phase structure becomes 
more complicated. Here we find that when the charge \q\ is greater than a critical value 
\q c \, there exists a globally stable black brane solution at every temperature in between 
zero and infinity. However, when \q\ < \q c \, there is a range of temperature, where there 
exist three black brane solutions with different radii. The black branes of the smallest 
and the largest sizes are locally stable as they correspond to the local minima of the free 
energy. On the other hand, the intermediate size black brane is unstable, corresponding to 
the maximum of the free energy, and never exists. The free energies of the largest and the 
smallest black branes are not the same and one is greater than the other which depends on 
the temperature of the system. However, there exists a transition temperature, for a given 
charge \q\ < \q c \, where the two free energies become the same. This is the temperature 
where two black brane states the smallest and the largest can coexist and can make a 
transition freely from one phase to the other just like the van der Waals-Maxwell liquid- 
gas phase transition. Since this transition involves an entropy change, therefore, it is a 
first order phase transition. This structure was also noticed in AdS [6, 7], asymptotically 
flat and dS black holes [9, 8] in canonical ensemble. When the temperature, in the range 

4 One expects this to be a topological first order phase transition since this involves a topology change 
as well as an entropy change. 

5 When the temperature is of the order of Planck scale, classical gravity description will break down 
and one must consider quantum gravity. 
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we mentioned, is less (more) than the transition temperature, the smallest black brane 
has lower (higher) free energy than the largest black brane. So, in one case the smallest 
black brane is globally stable and in other the largest black brane is globally stable. There 
is a first order phase transition line which depends on the charge \q\ < \q c \ and the line 
gets shrunk as we increase the charge ending up at a second order phase transition point 
(critical point) when \q\ = \q c \. We calculate the critical exponents at this critical point 
and found that they have very similar structure as in the black hole cases with a universal 
critical exponent [8, 6] for the specific heat as —2/3. Since the stability analysis for the 
black branes, in general, is quite complicated and extracting exact values of the parameters 
may not always be easy, we will illustrate the behavior by numerical calculations except 
for some special cases specified later. 

This paper is organized as follows. In section 2, we give the non-dilatonic black p- 
brane solution in space-time dimensions D and derive the action from which we study 
the stability and analyze the phase structure. The details of the general stability analysis 
is discussed in section 3. The phase structure for the case of zero charge in the cavity 
is considered in section 4. Section 5 discusses the case of non-zero charge. The critical 
exponents are given in section 6 and we conclude in section 7. The dilatonic brane cases 
are also mentioned and discussed in appropriate places in the respective section and how 
the whole analysis of the paper carries over to the dilatonic branes is discussed in the 
Appendix. 

2 Black £>-brane solution and the action 

The black p-brane solution was originally constructed as a solution to the ten dimensional 
supergravity containing a metric, a dilaton and a (p + l)-form gauge field [10]. This 
was generalized to arbitrary dimensions in [11]. These solutions are given in Lorentzian 
signature, but for the purpose of studying the thermodynamics [13] we will write the 
black p-brane solution (without the dilaton field as we will be studying the non-dilatonic 
branes, the dilatonic branes will be considered in the Appendix) in Euclidean signature 
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as (see for example [22]), 

ds 2 = A + Aj^~ 2 dt 2 + A^ 5Z(^) 2 + (A+A.)- 1 dp 2 + p 2 dQ\ +1 

dt A dx 1 • • ■ A dx p 



d d - X 



i=l 



d/2 / \ d/2' 



r+J V P 2 



~(r_r + ) 
P 

In the above we have defined, 



d/2 



F[ p+2 ] = <^4[p+i] = — *d — dp Adt A (ir 1 ■ • • A <ix p (1) 



— (7) « 2 » 

where r± are the two parameters characterizing the solution which are related to the mass 
and the charge of the black brane. The metric in (1) has an isometry S 1 x SO(d — 1) x 
SO(d + 2) and therefore represents a (d — 1) = p-brane in Euclidean signature. The total 
space-time dimension is D = d + d + 2, where the space transverse to the p-brane has the 
dimensionality d + 2. A p-brane couples to the (p + l)-form gauge field whose form and 
its field-strength are given in (1). It is clear from the Lorentzian form of the above metric 
that when r_ = 0, it reduces to D-dimesional Schwarzschild solution which has an event 
horizon at p = r + , whereas, at p = r_, there is a curvature singularity. So, the metric in 
(1) represents a black p-brane only for r + > r_, with r + = r_ being its extremal limit. 
Note that in the above we have defined the gauge potential with a constant shift, following 
[19], in such a way that it vanishes on the horizon so that it is well-defined on the local 
inertial frame. For the metric in (1) to be well defined without a conical singularity at 
p = r + , the Euclidean time in the metric is periodic with a periodicity 

which is the inverse of temperature at p = 00. The local 

P = Ay'Al^'/r, (4) 

which is the inverse of local temperature at p when in thermal contact with environment 
at the same temperature. For the canonical ensemble we have fixed p at the wall of cavity 

d(d-l) 

denoted by ps, fixed local temperature at ps, fixed local brane volume V p = At (D ~ 2) V*, 
with V* = f d p x and fixed charge defined as 



p 

V2kJ V2k 



Qd = -jU / *F [p+2] = ^d{r + r_.f' 2 . (5) 
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In (5) k = a/SttGd, where Gd is the D- dimensional Newton's constant, *F[ p+ 2] denotes 
the Hodge dual for the (p + 2) form- field given in (1). Also Q n denotes the volume of a 
unit n-sphere. With these data we will evaluate the action. 

The relevant action for the gravity coupled to a (p + l)-form gauge field in a manifold 
M of dimension D with the Euclidean signature has the form, 

Ie = i E (g) + Ie(F) (6) 
where the first term is the purely gravitational action, 

Ib{9) = ~ I d D x^R+± [ d D - l x^{K-K Q ) (7) 

lK JM K JdM 

consisting of the usual Einstein-Hilbert term and the Gibbons-Hawking boundary term 
[23] where DM denotes the boundary of the manifold M. In the above K is the trace of 
the extrinsic curvature K^ v defined as, 

= -- (V M n„ + V^) , and so, K = -V M n M (8) 

where is a space-like vector normal to the boundary and is normalized as n M n M = 1. 
Also 7 Q/ 3 is the boundary metric with a, ft the indices of the boundary coordinates and 
7 is its determinant. K is the subtraction term which serves as an infra-red regulator so 
that a finite result can be obtained [1, 13, 19, 3] 6 . This is calculated by embedding the 
same surface in the flat space. 

The second term is due to the form-field F[ p+2 ] and its expression in the canonical 
ensemble has the form, 

Note that for canonical ensemble the charge in the cavity is fixed. For a given gravitational 
configuration in Euclidean signature, the corresponding thermodynamical partition func- 
tion in the saddle point (or the zero- loop) approximation can be obtained as Z e~ lE , 

6 In the charged case, though this flat space is a solution of the corresponding equations of motion, it 
is not a thermal state and serves only as an infra-red regulator. One way to see the rationale behind is 
to first calculate the corresponding action in the grand canonical ensemble for which the flat space is a 
good reference state subtractor. Then the Euclidean action in the canonical ensemble can be obtained 
from the calculated Euclidean action in the grand canonical ensemble by a thermodynamical Lcgcndre 
transformation, noting that the respective action is related to /3 times the corresponding ensemble poten- 
tial to leading order. The canonical ensemble action so obtained is also finite and is actually equivalent 
to that by the subtraction procedure mentioned in the text. 
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where Ie is evaluated for the given configuration (see for example [20]) and we will eval- 
uate it in the black p-brane configuration given in (1). On the other hand, we also have 
Z = e~ l3F , where /3 is the inverse temperature of the ensemble and F is the corresponding 
free energy (for canonical ensemble this is Helmholtz free energy). Therefore, we have in 
this approximation F — Ie/P- This is the relevant quantity we are going to evaluate for 
studying the stability of the black branes and its phase structure. 

To evaluate Ie given in (6), we can use the equation of motion of the metric and 
obtain, 

d — d 



R 



2(D - 2)(d4 
Substituting (10) in the action we rewrite Ie as, 

d 



p 2 



(10) 



2n 2 (D-2)(d 
1 1 



H 



d u ~ Y x^{K - K Q ) 



2k 2 d\ 



d D - l x^n^F 



dM 



flfMX—Hd 



A 



Ml— Mci 



In the above we have used d + d = D — 2 as mentioned earlier and dd = 2{D — 2) for the 
non-dilatonic branes which is a solution of supergravity with maximal supersymmetry, 
the type we are considering. We will evaluate each term in (11) separately. For that 
purpose we need the forms of the normal vector n M , the trace of the extrinsic curvatures 
K and K which can be calculated from the metric in (1) and are given as, 



K 



(A+A^) 1 / 2 ^ 

V9 



d+1 
P 

d+1 
P 



(A + A_ 



,1/2 



d 



2p 



d+1 



a! 



1/2 



1/2 



(12) 



Now from the metric and the form-field given in (1), it is straightforward to evaluate the 
action (11) as, 



'^2^ Pb 



1/2 



A_ 



2 hr + d he +d(A + A 



A_ 



1/2 



.1/2 



2(d+ 1) 



(13) 



P=PB 



Note in the above that p is fixed on the cavity as p#. Now since we have Ie = /3F, with 
F = E — TS, the Helmholtz free energy in the canonical ensemble, so this implies that 
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Ie = (3E — S, with E the energy in the cavity enclosed and S, the entropy for the system. 
Following Braden et. al. [19] for black hole, we expect that the entropy in the present 
case (which is just that of the black brane under consideration) should be independent of 
the choice of the location of the cavity. However, this is not manifest in (13). For this we 
can rewrite (13) as, 



pv p n 



2k 2 
4vr VI* 



(d + 2) 



A \ 1/2 



2(d+l) 



2k 1 



P=PB 



(14) 



where we have used the expression of /3* given in (3). The last term in (14) is precisely 
the entropy of the non-dilatonic black p-brane 



S 



2k 2 + 



3+3 



2k 2 



1 _ 1 

2 d 



1 - 



(15) 



where we have used V„ 



d(d-l) 
^ 2(13-2) 



V* and is independent of the location of the cavity. We 



can also read off the energy for the cavity as 



E{ P l 



2k 



_-2 



Pb 



{d + 2) 



A. 



1/2 



rf(A+A_ 



,1/2 



2(d 



(16) 



P=Pb 



which approaches the ADM mass at pb — > oo. In casting the Euclidean action in the form 
of (14), we actually treated the fixed /3 at the wall of cavity, the inverse of temperature 
of the environment, to be independent of r + for the time being. Then the internal energy 
E and the entropy S are functions of r + only when pB,V p and Qd are all fixed as in this 
ensemble. The stability of the black p-branes can now be discussed from the action (14) 
when the second line of (15) for entropy is employed to which we turn in the next section. 



3 Generalities and stability of black ^-branes 

We have mentioned that in the canonical ensemble the charge inside the cavity given in 
(5) is fixed. This implies that r_ is not an independent parameter and can be expressed 



S 



in terms of r + as, 



V2nQ d 

n d+1 d 



2/d 



1 



(OS 



*N2 



;i7) 



where we have defined Q* d = [(y/2KQ d )/(n i+1 d)]^ d . With these, A_ will be erven as, 



r+P 



We are now considering the canonical ensemble for which the quantities /3, V p , ps and 
Qd are all fixed (this implies that /?*, V* are not fixed), and therefore the only parameter 
which can vary is r + . This implies that the entropy will change with r + , therefore if we 
vary the action (14) with respect to r + and set that to zero (that is, at the stationary 
point), we must be able to recover (3 (since the temperature is a conjugate variable to 
the entropy) and this in turn will give a non-trivial check that the form of action (14) 
we have obtained is indeed correct. For this we will replace V* by V p using the relation, 



V„ = A 



1 j d 

2" l "2(D-2) 



2 



V* = Af. d V* where we have used dd = 2(D — 2) in the second equality 



and so, (14) takes the form 
Ie = 



f3V p Vt d+l j 
Pb 



2k 2 



4tt V P A. 



(rf + 2) 



A_ 
A~ 



1/2 



rf(A+A_ 



1 1/2 



2(d 



" 2 ~ d Q 



2k 2 



d+l r d+l 




(19) 



In the expression of A± here and below the variable p must be replaced by constant ps 
the size of the cavity. From dlE/dr + = 0, we obtain after some simplification, 



l3d-ATir + A l ' 2 Aj d I 1 



i_i 

d 2 




0. 



(20) 



Since the second factor in the l.h.s. of (20) is greater than zero, we must have 



-^A 



V- 



i _i 

d 2 



d 



A_ 



A 1 ' 2 A_ 2(D " 2) (3* . 



(21) 



This is precisely the correct form of /3 given in (4), which we obtain from the metric in 
(1). We will use the relation (21) to discuss the stability and the phase structure of black 
p-branes. For this purpose let us rewrite /3 explicitly as, 

1/2 



P 



47rr_| 



i 



Q 



-.2d 



1 _ 1 

d 2 



,2(/ 



P d B 



ofX 

r%p d B ) 



(22) 
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where we have used (17). Now following refs. [19, 8, 9], we define 



x 



r+ 

PB 



< 1, 



P 



B 



PB 



(23) 



where the dimensionless parameters b, q are fixed 7 but x is the only parameter which 
can change. Note that the parameter b is related to the inverse of temperature of the 
environment, q is related to the charge and x is related to the horizon size. Also, as 
{Qd) 2 / r + = r ~/ r + < 1) x > Q- in terms of these parameters the above equation of state 
(or thermal equilibrium condition) for /3, (22), can be rewritten as, 



where 



b q (x) = - 
a 



b = b q (x) 
x i/d( 1 _ x y/2 

r 

i 



1 - si 



1 _ 1 

2 2 



1 - 21 



(24) 



(25) 



In the above (25) we have used 

A + = 1 - x, A_ = 1 - 
If we define the reduced Euclidean action as 



1 - 



X 



X* 



(26) 



2n 2 h 



-b 



—x 



1/2 



(d +2) 



x 



1 - 



+ rf(l-x) 1/2 ( 1 



X 



1/2 



- 2(d 




one can show 



dx 



f q (x)[b-b q (x)] 



(27) 



(2? 



where is as given above in (25) and 



-1/2 



1 - 



X 



'1/2 



d + 2 - 



d+2 1 



1 



> 0. (29) 



7 Since the charge as well as Q* d denotes the corresponding absolute value, so q also denotes the 
absolute value. Therefore, from now on, we use q to denote the absolute value of the reduced charge. 
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The equation of state (24) and the reduced action (27) are the only relevant equations 
which we need for the discussion of stability and phase structure for the black branes. 
While we derive them for the non-dilatonic branes, they remain true even for dilatonic 
branes as we will show in the Appendix. In other words, the analysis given in the following 
sections works for both the non-dilatonic and dilatonic branes. These two equations 
depend only on the parameters b, q, d and the variable x. The only relevant dimensionality 
is d which remains unchanged under the so-called double- dimensional reductions [21]. In 
other words, the black branes related by this reduction will have the same stability and 
phase structure at least in the approximation employed. For example, the D = 11 M2 
branes and the D = 10 fundamental strings share the same above mentioned properties 
and so do the other branes related by this kind of reductions (see, for example, [22]). 
Now from (28) we find that (at the stationary point) 

ox ox 1 
db ^<0 ^>0, (30) 



dx dx 2 

dx 



Therefore, the system will be stable (at least locally) when ^ < 0, corresponding to the 



minimum of the Euclidean action or the free energy. On the other hand, if ^ > 0, the 
corresponding Euclidean action is maximum and the free energy is also maximum and 
the system will be unstable. Since the phase structure for the black p-branes are quite 
different for the chargeless case and charged case, we will discuss them separately in the 
next two sections. 



4 Chargeless case 



In this section we will discuss the stability and the phase structure for the chargeless 
black p-brane. So, we will put q = in all the expressions we obtained in section 3. The 
expression for the reduced action now takes the form, 

I E = -2(d + l)b [(l-x) 1/2 - 



i+i 

X <*, 



whence we have 

OIe 

dx 

Where bo(x) and fo(x) are given as, 



fo(x) [b -fool 



X 



(31) 
(32) 



bo(x) = 4x^(1 -xf 2 , 
d 

f {x) = (d+i)(i-xr i/2 >o. 



(33) 
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Figure 1: The typical behavior of bo(x) vs x in the chargeless case. 



The equation of state is 

b = bo(x). 

At the stationary point of the action, we have 

d 2 I E 



dx 2 



dx 



(34) 



(35) 



where x is a solution of the equation (34). Since b (x) = at x — 0, 1 and b (x) > 
between < x < 1, so bo(x) has a maximum in between which can be determined from 
dbo(x)/dx = and has the value, 



T 



2d / d + 2 



47TPB 



at 



l/d 



PB 



(36) 



(37) 



c/ + 2 + 
So, there exists a maximum 6 max (or minimum temperature T min ) above (or below) which 
the system can not be in a black brane phase. Now since OIe/Qx > (from (32)) and 
Ie = at x = (from (31)), therefore the system favors the 'hot flat space'. Note that 
for the four dimensional black holes D = 4, and d = d = 1, and so, T m i n = 3v / 3/(87tpb) 
and r +max = (2/3)p£ match exactly with the values found in [13, 9]. We will have the 
same T m ; n and r +max for black strings in D = 5, black membranes in D = 6, up to black 
D6 branes in D = 10 since these branes are related to the four dimensional black hole via 
the double-dimensional reductions. 
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For b smaller than 6 max , we can have two solutions from the equation of state (34), 
but only the large solution x 2 (> x max ) will be locally stable given the relation (35). This 
behavior of bo(x) versus x is depicted in Figure 1. However, this does not necessarily 
imply that the system is in the black brane phase. Only when the local stability becomes 
a global one, then the system is indeed in the black brane phase. The corresponding x 2 
can be determined from requiring the action at the stationary point be negative. The 
action (31) can be expressed using (34) as, 

where we have defined 

y = y^]— (39) 

with X2 = x. So, the necessary condition for the global stability can be seen from (38) to 
be 

V < ~t~:- (40) 



This gives 



Now we find 



d + 2 



4(d+l) 

X > Xg (J_|_ Xnia - X ' \ I 




T^^iH^l) . (42) 

So only when < b < b g (in this case x g < x 2 < 1), the system is in the black brane phase. 
On the other hand, for b g <b < 6 max (in this case x max < x 2 < x g ), the system, though 
locally stable, will eventually tunnel to the 'hot flat space' at the same temperature. 
However, at b = b g , Ie = and so, both the black brane phase with x = x g and the 
'hot flat space' phase are possible. In other words, this is the place at which the two 
phases can coexist and the corresponding temperature is the phase transition one. This 
phase transition is both a topological and a first order one since both the topology and 
the entropy of the two phases are different before and after the phase transition. Note 
that when x — > 1, b — > as required from (34), but now b/y — > 1/d, so the action is 
still finite and is Ie = — 1, implying that the system is stable. For the four dimensional 
black hole we find that the temperature (T g ) where the large black brane becomes globally 
stable is 27/(327rps) which matches with the value found in [13, 9]. Also it is clear from 
Figure 1 that, as bo(x) decreases or the temperature increases the size of the small black 
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brane decreases and that of the large black brane increases and eventually when bo — > 
or T — > oo, the size of the small black brane goes to zero and the size of the large black 
brane approaches the size of the cavity. The phase transition we found in the present case 
is analogous to the Hawking-Page transition for the AdS black holes. 

5 Charged case 

In this section we will study the stability and phase structure of black p-brane in the 
more general case where the charge enclosed by the cavity is non-zero and fixed. In 
section 3, while discussing the generalities for the stability of black p-brane we found that 
the system will be stable when db q {x)/dx < and will be unstable when db q (x)/dx > 0. 
We will show that there exists a critical charge q c , such that when q > q c , the system will 
be globally stable as db g (x)/dx is always less than zero (see Figure 3 below), but when 
q < q c , db q (x)/dx > in some region as shown in Figure 2. We mentioned in section 3 
that the parameter x lies between 1 and q. Since b q (x) — > oo as x — > q and b q (x) — > 
as x — > 1 (see eq.(25)), therefore, when q < q c , b q (x) does not decrease monotonically (as 
seen from Figure 2) and db q (x)/dx > in some region of q < x < 1. In other words, 
we should have a minimum of b q (x) (fr mm ) occurring at x — x m [ n and a maximum of 
b q {x) (6 max ) occurring at x = x max . When b min < b q (x) < b 

max and X m [ n <C x <C x max , 

db q (x)/dx > and the system is unstable in this region and there exist no stable brane 
phases. On the other hand, when b g (x) > 6 max or b q (x) < b min , db q (x)/dx < and the 
system is stable (at least locally). So, for a given b with b min <b< 6 max , there will be 
three solutions to the equation of state (24). If we denote the three solutions as X\) x<2 
and x-$ with x\ < x<i < x%, then x\ and £3 correspond to the local minima of the free 
energy and X2 corresponds to the maximum. Among the two minima, one expects that 
the one with the lower free energy will be globally stable and a transition will occur from 
the state of higher free energy to the lower free energy. So, it is important to find which 
of the two black branes have lower free energy. To determine this we write from (28), 



M^3) - Ie(xi) = 5(x 2 ,xi) - S , (x 3 ,x 2 ), 



(43) 



where 




(44) 



(45) 
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Figure 2: The typical behavior of b q (x) vs x when there is a phase transition (q < q c ). 

with b q (x) as given in (25). Note that for b = 6 mm , the points £2 and x\ coincide and so 
S(x2, xi) = and S(x2, £3) takes a maximum value. Similarly, for b = 6 max , the points X2 
and X3 coincide and therefore, £(£2,23) = and S(x2,xi) takes a maximum value. Thus 
both the function S(x2,xi) and S(xs,X2) change continuously from to their maximum 
value as b is varied from b min to b max . So, in between there must exist a fixed value of b 
which we call b t (inverse of which is related to a phase transition temperature) for each 
given charge q < q c , for which S(x2,xi) = S(xs,X2) and therefore the Euclidean action 
or the free energies of the two stable black brane configurations of sizes x\ and £3 are 
the same. In other words, this is a phase transition temperature where the two black 
brane phases coexist and make a transition freely from one phase to the other, much like 
a van der Waals- Maxwell liquid-gas phase transition. This was also noticed earlier for 
asymptotically AdS, dS and flat black holes in canonical ensemble [6, 9, 8]. 

Now since at b max , S(x 3 ,x 2 ) = and S(x 2 ,Xi) is maximum and at b t , they are the 
same, so, if b lies in between i.e., 6 max >b>b t , S(x 2 ,Xi) > S(x 3 ,x 2 ), or Ie(x 3 ) > Ie{xi). 
So, in this case the small black brane phase is globally stable. Similarly, since at 6 m i n , 
5'(a;2,^i) = and £(£3, £2) is maximum and at bt, they are the same, so, if b lies in 
between i.e., 5 mm < b < b t , S(x3,X2) > S(x2,xi), or Ie(xi) > Ie(xs). So, in this case 
large black brane phase is globally stable. Thus we conclude that if the temperature is 
above the transition value, but below T max , large black brane phase is globally stable and 
if it is below the transition value, but above T m j n , the small black brane phase is globally 
stable. So, there will be a phase transition from the small black brane to large black brane 
or vice versa depending on whether the temperature of the black brane is below or above 
the transition temperature. 
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Figure 3: b q {x) decreases monotonically with x when q < x < 1 (q > q c ). 

It is clear from the expression of entropy given in the last term of (14) that for given 
q < q c , the entropy will depend on the parameter r + or x and so, the entropy will change 
during the phase transition we just mentioned. Therefore, this is a first order phase 
transition for which the entropy has a discontinuity. In fact there is a first order phase 
transition line when we move the charge q < q c towards q = q c and the line gets shrunk 
ending up at a second order phase transition point (critical point) for q = q c , which occurs 
at x = x c , and where the entropy discontinuity disappears. 

Having understood the qualitative features of the equilibria and the phase structure 
of black p-branes with non-zero charge inside the cavity, we will try to understand the 
structure in a more quantitative way and then corroborate our observations by numerical 
calculations how the various situations for q > q c , q = q c and q < q c described here arise in 
three different values of d, for examples, for the non-dilatonic branes (we will consider only 
d = 3 in D = 11, which corresponds to M5-brane, d = 6 in D = 11, which corresponds to 
M2-brane and d = 4 in D = 10, which corresponds to D3-brane) though these calculations 
work for the corresponding dilatonic branes related via the double-dimensional reductions 
in other dimensions as well. 

For understanding the stability, as we mentioned, the quantity to look at is db q (x)/dx. 
From (25) we find, 



96 (a.) x 1 '* {(1 + l)x 4 - [1 + (2 + i)q 2 }x 3 - 3g 2 (| - l)x 2 + q 2 [d - 1 + f q 2 ]x - d ? 4 } 



dx 




(46) 
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The position of the extremality will be determined from the vanishing of the numerator 
of the above equation (46), i.e., 





x 




x 2 + q 2 d-l + —q 2 x-dq A = 0, 

(47) 

This is a quartic equation and has four roots in general. We will make some observation 
about the roots of this equation which will support the various structures we described 
qualitatively in this section and then give some numerical solution of this equation in 
some special cases. First note that the discriminant of the above equation (47) which 
tells us about the roots has the form, 



AM) 



(<? 2 



i)Y 
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(4(d - 1) - 3d(4 + d)q 2 f - 108rf 2 (2 + d-d 2 ) 2 q 2 (l - q 2 ) (48 



The discriminant will vanish within < q 2 < 1, if the last factor within the square bracket 
in (48) vanishes in that range. This will be determined by the intersection point of the 
first term, i.e., a cubic curve with the second term, i.e., a parabola. The parabola meets 
the g 2 -axis at q 2 = and at q 2 = 1 and remains positive in this range. On the other hand 
the cubic curve takes a positive value of 4 3 (d — l) 3 for d > 1 at q 2 = 0, monotonically 
decreases and meets the g 2 -axis at q 2 = 4(d — l)/(3d(4 + d)) < 1 (for d > l) 8 . Therefore, 
there is a unique crossing point q 2 of the cubic curve and the parabola in the range 
< q 2 < 1 with < q 2 c < q 2 = 4(d - l)/(3d(4 + d)) < 1. This shows the existence of a 
unique critical point at q — q c , where db q (x)/dx vanishes. Note that since this is a single 
extremum, as b q (x) varies from oo to 0, this can not be a maximum or minimum, but is 
an inflection point where d 2 b q (x)/dx 2 also vanishes. This feature is reflected in Figure 4 
below. Note that for d = 1, q 2 = and the intersection now occurs at go — which is 
not in the range of < q 2 < 1, therefore we don't expect a critical point and further a 
possible phase transition to occur. This will be checked explicitly in a subsection of this 
section later. 

When q 2 > q 2 , we have A(q, d) < 0, so (47) (since this is a quartic equation, there are 
four roots of this equation in general), must have a pair of complex conjugate roots. Also 
since the ratio of the coefficient of x 4 term and the constant term (it is — dq 4 /(l + d/2)) 
is negative, implying that the product of the four roots are negative, so one of the roots 
must be negative. Hence, there can be at most one root in the region q < x < 1. However, 
for d > 2, as the l.h.s. of equation (47) is positive at both x — q (it is (— 2 + d)(— 1 + q 2 )q 3 ) 
and x = 1 (it is d(—l + q 2 ) 2 /2), the number of roots, if existing at all in the region 



3 We will discuss d — 1 case separately in a subsection of this section. 
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q < x < 1, must be even 9 . So, there can not be any root of (47) in the range q < x < 1 if 
q 2 > q 2 c and this is consistent with the typical behavior given in Figure 3. 

When q 2 < q 2 , we have A(g, d) > and this implies that (47) has either two pairs of 
complex conjugate roots or all real roots. Also since the product of all four roots must be 
negative, then the roots must all be real and the number of negative roots must be odd 
(one or three). By the argument of previous paragraph, the number of roots, if existing at 
all in the range q < x < 1, must be even (two or none). The finding that x c falling in the 
region q < x < 1 for q = q c , for each of the cases considered in the following, implies that 
indeed there exist two roots in the region q < x < 1 for q 2 < q 2 . This is also consistent 
with the typical behavior given in Figure 2. 

Now we give some numerical calculation to illustrate the above picture. When A(g, d) = 
0, one expects that q 2 in (48) has only one real positive root q 2 , For this critical q c , we 
should have x m i n = x max in the region q c < x < 1. Let us examine carefully if this is 
indeed true. Let us first consider M5-brane, i.e., d = 3 first. Now A(q c ,d = 3) = gives 
us, 

-512 + 2764% - 11080% 2 + 250047?/ 3 = 0, (49) 

where we have defined y = q 2 . One can check that indeed y has only one positive real 
solution given by 

4 (l!4 - 4478 ( _ 450413+ 2 222607 ^ ) 1/3 + 2^ (-450413 + 222607^) 1/3 



y 



3087 

0.020058, (50) 



whence we get 

q c = 0.141626. (51) 

Substituting this in (47) we obtain 

-0.00120697 + 0.0419264x - 0.030087x 2 - 1.0702x 3 + ^x 4 = 0, (52) 

which indeed gives a unique solution in the region q c < x < 1 as, 

x c = 0.292675 (53) 

(the other two irrelevant solutions are x = —0.187352 and x = 0.030083). So, we now 
have the critical b c as, 

xl /3 (i-x c y/ 2 , . 

0.199253. (54) 



3(1 - |)V6(1 - |)V3 
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Figure 4: The typical behavior of b q (x) vs x when there is a turning point for which 

•^min -^max %c (.9 9c)- 

From (47), one expects that there exist a minimum of b q (x), occurring at x — x m i n and 
a maximum, occurring at x — x max , when q < q c , corresponding to Figure 2 (A(q, d) > in 
this case), while such x min and x max should not exist when q > q c corresponding to Figure 
3 (A(q r , d) < in this case). Let us take two explicit examples, one for each case, showing 
that this is indeed true. When we take q = 0.150000 > q c which is slightly larger than 
the critical value, we find indeed that there exist no real solution in the region q < x < 1. 
While if we take q = 0.130000 < q c which is slightly smaller than the critical value, we 
find that indeed there exist two solutions in the region q < x < 1, one is x m \ n = 0.229179 
and the other is x max = 0.341762, as expected. 

For M2-brane, d = 6, we find exactly the same behavior. In this case A(q c , d = 6) = 
gives (see (48)) 

-8000 + 3264192y - 4992192y 2 + 5832000|/ 3 = 0, (55) 
where again we have defined y = q 2 c . The only real positive solution of (55) gives, 

321 " 9 5 06 ( - 2 03i 67+ 7 i 6825 o^ ) 1/3 + 272/3 (-203167+ 168250 V2) V3 

V = 

y 1125 

w 0.00246007, 
which gives 

q c = 0.049599. 
9 We will discuss d = 2 in a separate subsection of this section 



(56) 
(57) 
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This in turn gives the unique solution 



x c = 0.175176, (58) 
from the corresponding equation (47) 

-0.0000363115 + 0.0123548x - 0.0147604x 2 - 1.0123x 3 + 4x 4 = 0, (59) 
in the region of q c < x < 1. We also have 

b c = 0.116698. (60) 

The other cases q > q c and q < q c can be discussed similarly as above. 
For D3-brane, d = 4, A(q c , d = 4) = gives 

-1728 + 214272y - 504576?/ 2 + 884736y 3 = 0, (61) 

whose unique positive real solution is 

1/3 , .2/3/ o„o^ , /FX V3 



y 



73 ~ 1315 ( Wie^yls ) + 52/3 (~ 34367 + 16512 v / 6) : 

384 

0.00822139, (62) 



which gives 

q c = 0.090672. (63) 
With this critical q c , we have the unique 

x c = 0.238800 (64) 

in the region of q c < x < 1 from the corresponding equation (47) 

-0.000270365 + 0.0250697x - 0.0246642a; 2 - 1.03289x 3 + 3x 4 = 0. (65) 

Now 

b c = 0.159921. (66) 

Once again, the other cases q > q c and q < q c can be similarly discussed. 

We have given the numerical results for d = 3, 6, and 4 since they are related to 
the M5, M2 and D3 branes. As indicated earlier and with the discussion given in the 
Appendix for dilatonic branes, any brane, non-dilatonic or dilatonic, related to each of 
these branes by the so-called double-dimensional reductions will share the same properties, 
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since the only dimensionality entering this discussion is d and it remains the same under 
this reduction. 

Note that since we have 1 < d < 7, for completeness and to show that the critical value 
x c always falls in the range q c < x c < 1, we give the results also for d = 5 and 7 (d — 1, 2 
will be discussed separately). The critical values for d — 5 are (q c = 0.064944, x c = 
0.202012, b c = 0.134632) and for d = 7 they are (q c = 0.039529, x c = 0.154691, b c = 
0.103210). We collect the relevant quantities for the cases 3 < d < 7 in the tabular form 
as, 



d 


Qc 


X c 


b c 


3 


0.141626 


0.292656 


0.199253 


4 


0.090672 


0.238800 


0.159921 


5 


0.064944 


0.202012 


0.134632 


6 


0.049599 


0.175176 


0.116698 


7 


0.039529 


0.154691 


0.103210 



The table above shows that the critical size of the black brane always lies in the range 
q c < x c < 1 , where q c is related to the absolute value of critical charge and b c is related to 
the inverse critical temperature as defined earlier. We observe that the critical quantities 
all decrease as d increases. 



5.1 d = 2 

The general equilibria and the phase structure that we discussed in this section does not 
apply to d — 2, 1 cases as we mentioned before and so, we will study these two cases 
separately here in the next two subsections. Let us first discuss d = 2 case. We find from 
(25) that for d = 2 

, , , x^il-x) 1 ' 2 

M x ) = — ~ 77TJ2- ( 67 ) 



The corresponding equation (from (47)) for finding the extrema for d = 2 is 

2x 4 - (1 + 3g 2 ) x 3 + q 2 (l + 3q 2 ) x - 2q 4 = (68) 
which can be factorized as, 

{x 2 - q 2 ){x - x+)(x -x-) = (69) 
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Figure 5: The typical behaviors of b q (x) vs x for d = 2. 

where, 

x± = \(i + 3q 2 ± V^) (70) 

with A = (1 — q 2 )(l — 9q 2 ). The x 2 = q 2 solutions are irrelevant here and x± can be real 
only if A > implying < q < 1/3. Note that the discriminant for the present case from 
(48) has the form, 

A(g,2) = 4g 6 (l-g 2 ) 3 (l-9g 2 ) 3 (71) 

This gives the same requirement as A for the reality of the solutions. Comparing with 
our earlier discussion of A(q,d), we note that here there is no parabola and the cubic 
curve meets q 2 axis at go = 1/3- So, this is also the critical point q c = qo = 1/3 which 
is consistent with the relevant solution of A = 0. Further we note that for this case, 
b q (x = q) = \/g/2 which is different from oo (for non-zero charge), the value for d > 2, 
but the other end b(l) = remains the same. 

Now let us see the phase structure in detail (see Figure 5). When 1 > q > q c = 1/3, 
we see that there exists no extrema for q < x < 1, since now A < and 

db g (x) _ bgjx) 
dx 2 

So, for now, when < b < b q (q), we have a stable black brane 10 . For b > b q (q), there is 
no stable black brane and we don't have a description available for such phase. 

For q = q c = 1/3, A vanishes and the two roots x± are equal and has the value 1/3 = q c 
. So, b q {x = x c ) = b c is a fake critical point which is not accessible since this is also an 
extremal point. This point is only marginally stable so long as the thermodynamics is 

10 Unless we have a phase like 'hot flat space', but now carrying a charge, whose free energy can be the 
lowest to be the globally stable phase, this black brane phase is globally stable. 



2x 2 - x(l + 3q 2 ) + 2q 2 
x(l — x)(x — q 2 ) 



< 0. 



(72) 
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concerned. Now as we increase x beyond this value, b q (x) monotonically decreases and 
goes to zero at x = 1 and in this range db q /dx < 0. Thus we find that in the range 
b c = l/(2y / 3) > b > b q (l) = 0, there exists a stable black brane phase. However, if 
b > 1/ (2\/3) there are no black brane phase or such description is not available . In terms 
of temperature, if the temperature T is below y/3/ (2itpb), there are no black brane phase, 
however, if the temperature is in the range a/3/ (2hpb) <T < oo, there is a stable black 
brane phase. 

For q < q c = 1/3, A > 0, so in this case we have two real solutions of db q {x)/dx = 
given by (70). One can check that x_ < q and q < x + < 1 and so only x + is the relevant 
solution which gives a maximum of b q (x) at 

(l + 3g 2 + V(l-? 2 )(l-9g 2 )) V2 - v^v) r q 
& max = ^ ^=-± L > b q (q) = Vl. (73) 

Note that for q < x < x + , db q (x)/dx > while for x + < x < 1, it is less than zero. 
Therefore in the range < b < b q {q) = y/q/2, there is only one black brane phase with 
1 > x > x + which is stable (we assume that there exists no other stable phase) and in 
the range b q (q) <b< 6 max , there are two black brane phases in which the smaller one is 
unstable and the larger one is stable. For b > 6 max there is no black brane phase or such 
a description is not available. 

5.2 d= 1 

For d — 1, we find from (25) the form of the parameter b q (x) as 

1 

„2\ 2 



b q (x) = — ^ '—. (74) 
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From b q (x = q) = b q (x = 1) = and the fact that b q (x) > for q < x < 1, there must 
exist one and only one extremum which corresponds to a maximum, denoted as 6 max (see 
Figure 6). This is due to the absence of a critical point as discussed earlier for the present 
case. Also as discussed in the previous section, since now q$ = 0, so 



A(g,l) 



(15) 3 (g 2 - l) 3 g 



3„8 



16 P 16 

q q H 

H 125 y 125 



< 0, (75) 
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for < q 2 < 1. So there must exist a pair of complex conjugate solutions for the following 
extremal equation of b q (x) which is a special case of (47) for d — 1, 

~x 4 - (l + \A x 3 + ~ gV + \q' x - g 4 = 0. (76) 
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Figure 6: The typical behavior of b q (x) vs x for d = 1. 

Furthermore, as the product of four roots of the above equation is less than zero, this 
implies that the other two roots must be one negative and one positive. Given that there 
exists a maximum, this positive root must be the one we expect and should fall in the 
region q < x max < 1. Let us take a special case of q = 0.50 as an example. Now the four 
solutions of (76) are 

xi = -0.30, x 2 = 0.33 - 0.29a, x 3 = 0.33 + 0.29i, x 4 = 0.73, (77) 

which are exactly as expected with only x 4 = 0.73 > q as the location of maximum of 
b q (x) for this q = 0.50. The situation here is similar to the chargeless case but we have 
here q < x < 1 for b q (x) vs x. In addition, we do not have the corresponding transition to 
'hot flat space', i.e., the Hawking-Page transition, since the charge is fixed for the present 
system. 

For a given < b < 6 max , there are two blackbrane solutions: the small black brane 
is unstable while the large black brane is at least locally stable. Unless some new phase 
carrying the same charge is found with lower free energy, we can not justify whether such 
a large brane is globally stable or locally stable. 

This therefore completes our analysis of the equilibria and the phase structure of 
black p-branes in the case of non-zero charge. We have seen that the d = 2 serves as 
a borderline which distinguishes the d = 1 case from the d > 2 cases. For the former 
case, the introduction of charge doesn't change qualitatively the stability behavior except 
for putting a new lower bound on the x set by the charge. In addition, there doesn't 
appear to exist the obvious analog of Hawking-Page type transition. Apart from this, the 
phase structure remains basically the same as the chargeless case. For the latter case, 
however, the introduction of charge does significantly change the stability as well as the 
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phase structure from the chargeless described in detail in this section. One of 

striking features is the appearance of a critical charge which determines both the stability 
behavior and the phase structure of the underlying system. When the charge is less than 
the critical charge, there exists a first order phase transition line which ends at a second 
order phase transition point which is the critical point at the critical charge. We would 
like to remark that although the details of the phase structure is quite different for d > 2 
for black branes, the qualitative structure is very similar to those of asymptotically AdS, 
dS and flat black holes in canonical ensemble studied earlier [6, 8, 9]. 



6 Critical exponents 

We observed in section 5 for the case of charged black branes, that when the charge is 
below certain critical value, q < q c , there exists two stable black brane states in certain 
range of temperature. Let us denote the sizes of the two black branes as x s and x#, where 
the former denotes the small black brane and the latter is the large black brane. There is 
a transition temperature where the free energies of the two black branes are the same and 
at this temperature these two black brane phases coexist. This transition temperature, 
which is completely determined by the given charge q < q c , therefore, forming a first order 
phase transition line, can be described by T t (q) (with the subscript 't' denoting it as a 
phase transition temperature). As the charge increases, this phase transition line ends 
up to a second order phase transition point (critical point) at q = q c . We would like to 
calculate the critical exponents at this critical point. Expanding b q (x) around the critical 
point x c we have, 



1 &b n 



+ 



(n 



3! dx 3 

Note that at the critical point the first and the second order derivatives of b q (x) with 
respect to x are zero. From (25) we find, 



d 3 b n 



dx 3 



i 

6 Xc 



1 + - 



2 + 1 



d 2 xi(l 



1-% 



d-2 
2d 



1- £L 



2+1 



(79) 



where we have used 



^bg 

dx 



0. 



dx 2 



0. 



(80) 



We mentioned in section 2, that the form of the entropy can be read off from the last 
term of the Euclidean action (14). If we now define the reduced entropy for fixed ps and 
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Vp as, S = 2K 2 S/{Vt^ +l V p p d B 1 ), then from (27) we read off the form of reduced entropy 



as, 



d+1 / Q 

S = Attx^t [ 1 - *- 



X 



1_1 

"2 d 



2 1 d 



We now expand the entropy around the critical point as, 



S-S a 



dS_ 

dx 



[X - X c ) + 



(82) 



where from (81), 



dS_ 

dx 



2nx~a- 



2(d + l)x 3 - (3d + 4)gV + 2q 2 x + q 4 d 



dx 3 1 - 21 



2 \ 2^ d 



1 - ^ 



d-2 
2 \ 2d 

XT* 



J3) 



We would like to have the expansion of entropy not around the x c , rather around the 
reduced critical temperature r c = l/b c . Note that near the critical temperature, 

1 1 



1 1 d 3 b n 



b 2 3! dx s 



[x - x c y + 



M) 



where we have used (78). So, using (84) we can write (82) in terms of the reduced 
temperature as, 



C _ Q — 

ax 





hi 




1 d 3 6, 




3? — Xc 


3! dx 3 


X — Xc 



1/3 



T - 7V. 



1 1/3 



55) 



The reduced specific heat therefore can be calculated as 

_ dS_ 

1 



3 <9x 



3! dx 3 



1/3 



(r - r c )~ 2 / 3 + 



?6) 



Therefore, the critical exponent a of is —2/3 for black branes when d > 2. 

Let us consider M5 as an example. We can expand (24) around the critical value as 

1 d 3 b n 



b q - b c 



[x - x c )° + ■ 



3! dx 3 , 4 
-1.88444 (x-x c ) 3 + 



57) 
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S-S, 



c 



12.6268(x -x c ) + - 



(88) 



We have now 



T -T, 



c 



47.465(x - x c ) 3 + ■ 



(89) 



S-S c 



3.48741 (r-r c ) 1/3 + 



(90) 



and 



c v = 5.83415(r - r< 




(91) 



Note that the critical exponent for the specific heat has a universal value —2/3 as was 
also noted for the asymptotically AdS, dS and flat black holes earlier in [6, 8]. 

As indicated above already, the present analysis works for both the non-dilatonic and 
the dilatonic branes. 



To conclude, in this paper we have studied in detail the equilibrium states and the phase 
structures of the asymptotically flat non-dilatonic and the dilatonic black branes in a 
cavity in arbitrary space-time dimensions D. Although we mostly concentrated on the 
non-dilatonic branes, the whole discussion applies also to the dilatonic branes as we give 
the details of how the analysis can be carried over in the Appendix. We considered 
only the canonical ensemble and so the charge inside the cavity and the temperature at 
the wall of the cavity were held fixed. We employed the Euclidean action formalism to 
compute the thermodynamics and the phase structure of the black branes. There is a 
marked difference in the phase structure when the charge enclosed in the cavity is zero 
and non-zero. When the charge is non-zero, there is also a qualitative difference in the 
phase structure when the d > 2, d = 2 and d < 2. We discussed them each separately. 

For zero charge we found an analog of Hawking-Page transition even for these higher 
dimensional black objects in asymptotically flat background. So, we found for the zero 
charge case that there exists a minimum temperature given in eq.(36), below which there 
is no black brane phase and here the system will be in 'hot flat space' phase. But above 
this temperature there exists two black brane phase with different radii. The smaller 
one is unstable and the larger one is locally stable. When the temperature lies between 
the minimum value (36) and the value given in (42), the locally stable black brane will 
eventually tunnel into 'hot flat space' since the latter configuration in this region has 
lower free energy. There is a phase transition temperature given by (42) at which the 
black brane of size x g given in (41) and the 'hot flat space' can coexist. But above this 
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temperature, the larger black brane becomes globally stable and the system can remain 
in this black brane phase. However, the smaller black brane still remains unstable. As 
we increase the temperature the size of the smaller black brane becomes smaller and that 
of the larger black brane becomes larger. As the temperature tends to infinity the size of 
the smaller black brane tends to zero and the larger black brane approaches the size of 
the cavity. 

When the charge enclosed in the cavity is fixed but non-zero, we found that there 
exists a critical charge q c at or above which there is a single globally stable black brane 
phase in the absence of an analog of 'hot flat space' or some other unknown configurations 
with favorable free energy for the charged system. When the charge is below this critical 
value, then there exists a certain range of temperature denoted by T min and T max below 
and above which there is a single globally stable black brane phase. But within this range 
there are three black brane phase with different radii. The largest and the smallest of 
which are locally stable as they correspond to the local minima of the free energy, but the 
intermediate one is unstable as it corresponds to the maximum of free energy. We found 
that the values of the free energies of the black brane phases depend on the temperature. 
In fact there exists a unique transition temperature at which the free energies of the largest 
and the smallest black brane become equal and so these two phases can coexist at this 
temperature and can make a phase transition freely from one phase to the other much like 
the van der Waals-Maxwell liquid gas phase transition. Above this transition temperature 
the larger black brane is globally stable and below this temperature the smaller black brane 
is globally stable. So, there is a first order phase transition from smaller black brane to 
larger black brane or vice versa above or below the transition temperature as the entropy 
of the system changes in this transition. In fact when the charge increases from q < q c 
towards q = q c , there is a first order phase transition line which eventually ends up in a 
second order phase transition point (critical point) at q — q c . At this critical point we 
have calculated the critical exponents and found that the critical exponent of specific heat 
has a universal value —2/3. We have elaborated this phase structure both analytically 
and numerically to illustrate the various situations. We found that this general phase 
structure is valid only for d > 2, where d is related to the dimensionality of the black 
p-brane by d = D — p — 3. d — 2, 1 case has been considered separately in subsections 
5.1 and 5.2. d — 1 case is very similar in structure to the zero charge case except here b q 
vs x curve starts from x = q instead of x = as in zero charge case. Also here there is 
no analog of 'hot flat space' since the system has non-zero charge. For d = 2, we found 
that there exists a critical charge q c above which there is a single globally stable black 
brane phase when the temperature of the system is above a certain value, but below this 
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value of the temperature there is no black brane phase and we do not have a suitable 
description for this phase. When the charge becomes the critical value q c , again there 
exists a temperature T = \f3/(2ixp B ) above which there is a single globally stable black 
brane phase and below this temperature there is no description available for the charged 
system. When the charge is below q c , there is a certain range of temperature where there 
are two black brane phases, the smaller one is unstable and the larger one is globally 
stable. Above this range, there is a single globally stable black brane phase and below 
the minimum temperature we do not have a description available. This whole analysis 
works for both the non-dilatonic and dilatonic branes. 

The only dimensionality which is relevant to the stability and phase structure is d 
and this implies that the branes related via the so-called double-dimensional reductions 
have the same stability and phase structure at least in the leading order approximation 
adopted. For example, this implies that the D = 11 M2 brane, the D = 10 fundamental 
string and the D = 9 0-brane all have the same stability and phase structure, so do the 
D — 11 M5, the D = 10 D4, the D = 9 3-brane, upto the D = 6 0-brane, and so on. We 
also observe that the critical quantities (q c , x c , b c ) all decrease when d increases from 2 to 
7. 
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Appendix 



Here we will consider the dilatonic black p-brane solutions in D space-time dimensions. 
For studying thermodynamics we give their form in the Euclidean signature as, 



- d-l 2 
a , 



ds 2 = A+A_ D - 2 dt 2 + A°- 2 J2( dxl ) 2 + A + lA - d dp 2 + p 2 Al d n 2 ^ 



d+V 

i=i 



4r _ 1 'p a< / > 0/2 



d/2 / \ d/2 

~7~ 



dt A dx 1 A ... A cfe p , 



Fu+2] = dAu+ij = -ie a(f)o/2 d ( r V ^ 1 dp A dt A dx 1 A ... A dx p , 

pd+i 

e 2(4>-M = A a ; (92) 

where A± are as defined before in section 2. is the dilaton and O is its asymptotic 
value and related to the string coupling as g s = e^°. a is the dilaton coupling given by 11 , 

fl2 = 4 -^2 (93) 
In the metric given in (92), the Euclidean time is periodic with periodicity (3* given as, 

This is the inverse temperature at p = oo. The local /3 is given as, 

= a^aZ^+^/T (95) 

which is the inverse of local temperature at p. The black p-brane will be placed in a cavity 
with its wall at p — ps- It is clear from the metric in (92) that the physical radius of the 
cavity is 

p B = Al d p B , (96) 

while pb is merely the coordinate radius. It is this ps which we should fix in the following 
discussion and not ps (as in the non-dilatonic case). Also we fix the dilaton on the 
boundary which is the requirement of obtaining its standard equation of motion from 
the action given later. In other words, we fix the dilaton at p~B, which indicates that 



n Note that the form of a 2 is fixed by supersymmetry, in the sense that these are solutions of super- 
gravity with maximal supersymmetry. 
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the asymptotic value of the dilaton is not fixed for the present consideration and this is 
crucial for our discussion. By this argument we also have 

f± = Al d r± (97) 

and r± are the proper parameters which we should use in the present context. In terms 
of the 'barred' variables A± remain the same as before, 

A± = 1 - -± = 1 - -4 (98) 

PB P B 

For the canonical ensemble we have fixed local temperature at the wall of the cavity, fixed 

rf(ri-l) 

local brane volume V p = A5°~ 2) V* and fixed charge defined as, 



Q d = J^J e -«» * F [p+2] = ^=e-^d(r + r. 



\d/2 



= 5^e-^ 2 (f + f„) J / 2 (99) 

where in the last line we have expressed the asymptotic value of the dilaton by the fixed 
dilaton = <P(pb) at the wall of the cavity from the relation given in (92) and then 
expressed r± by f± from (97). 

With these data we will now evaluate the action. The relevant action for the gravity 
coupled to the dilaton and a (p + l)-form gauge field with the Euclidean signature has 
the form 

i E = i E (g) + i E (<P) + i E (F) (ioo) 

where, Ie{.q) is the gravitational part of the action, Ie{4>) is the action for the dilaton and 
I E {F) is the action for the form-field and are given as, 

I E {g) = ~f d D x^R+\[ d D - l x^y(K-K ), 

ZK JM k JdM 

I d D - l x s fin„e-" d >*F""'«-'»A fllla ... l , d , (101) 

JdM 



1 1 



2K 2 d\ 

The various quantities in the above actions have already been defined in section 2. The 
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equations of motion following from the action (100) have the forms, 



□c 



V ('e -a ( d ^_F m "' Md+lN 



in 



-d^d^ - -(d(f)) 2 g^v 
1 1 / 

__i p -a{d)4> I rp p MlM2-"Atd 

2d\ I r m»-mu li r i' 

a ( d ) r -a{d)4> F 2 

2(d + l)! 





2RTT)^ F 7^ 02) 

(103) 
(104) 



Using the equation of motion, the action can be reduced to: 

d 



2(D-2)K 2 (d + 
1 1 



— / d D x^e~ a ^Fl l + -J d D ^x^(K - K ) 

L )- J M K JdM 



2k 2 d\ 

From the metric in (92) we have 



JdM 



1*2— • 



1/2 A 3 4d tf M 



The extrinsic curvature for the p-brane can be calculated as before, 



K = -V/ = --9,(^") 
i 1 ° 2 1 / o 2 d 



d a 2 



A 



-i 



The extrinsic curvature can be calculated as 

(J+l)Al ! 



d + 1 



(105) 



(106) 



(107) 



(108) 



where we have defined p = At d p and note that with this redefined p, Kq takes exactly 
the same form as in the non-dilatonic case. 

Now calculating each term in the action (105) separately as before we obtain, 



2k 2 



A \ V 2 /A \ 1/2 

2(=) +rf (=H +rf(A + A_) 1 / 2 _2(rf + l) 



(109) 



Note that in the above action everything is expressed in terms of the 'barred' parameters 
defined earlier instead of the original parameters. Comparing the reduced action (109) 
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with the corresponding reduced action for the non-dilatonic branes (13) we find that they 
have exactly the same form in terms of the redefined parameters. Once we have this form 
of the action (109), we can rewrite it as before in the form 1% = j3E — S as, 



WW,- 



2k? 



(d + 2) ( £±J + d~(A + A_) 1 / 2 -2(d~+l) 



4* ,+i, 

2k? r + 



. 1 + i 



P=Pb 



(110) 



where we have used 



0* 




We can thus identify the entropy 



S 



V d + l r d+l I j 



2 + d 



2k 



.•2 



111^ 



(112) 



and the energy of the cavity as, 



AN 1 / 2 

(d + 2) ( ^ij + d(A + A_) 1 / 2 - 2 (d+l) 



(113) 



Note that the entropy has exactly the same form as that of the non-dilatonic brane and 
we find that the energy approaches the ADM mass at pB — > oo as expected. We would 
like to remark that all the quantities like energy, entropy and temperature in the cavity all 
have the invariant forms in terms of 'unbarred' (non-dilatonic case) and 'barred' (dilatonic 
case) coordinates. Energy expression given in (113) has already the same form as can be 
compared with (16). Entropy given in (112) can be written as, 



S 



2k 2 



2k 2 




(114) 



Comparing (114) with (15) we find that indeed they have exactly the same form. Similarly 
we have from (95) and (94) the expression for inverse temperature as, 



Airr + . i . -i 
— +^AiA d 



(115) 



33 



Comparing (115) with (21) we again find that they have exactly the same form. This 
therefore indicates that the thermodynamical quantities of the non-dilatonic branes es- 
sentially have the same structure as the dilatonic branes with the new physical parameters. 
Now using the expression of charge (99) we can write f_ in terms of f+ as, 



v^ e aH 2 y i = m_ (116) 

Where we have defined Q* d = {(V2nQ d e a & 2 )/(tt (i+1 d)} 1 / d . Note that since the charge Qd 
is fixed inside the cavity and the dilaton <j> is fixed at the wall of the cavity Q* d is also 
fixed. Therefore, f_ is not an independent parameter, but is dependent on f + as given in 
(116). Now using (116) we can write (115) as, 

p = ^( 1 _9^) H (i-iY'(i-3t) i , air, 
d \ rf ) { pi) \ f%- 9 %) 

Now as before we can define 



In terms of these parameters (117) takes the form 



i i 



Va-xWi-^'fi— V^w*> (I") 



d V x J V x . 

Comparing with (24) and (25) we find that this equation (119) has exactly the same form 
as for the non-dilatonic branes. The expression of b q (x) was crucial for our analysis for the 
equilibria and stability structure of the black p-branes. Since the dilatonic branes have the 
same expression for b q (x) as the non-dilatonic branes, the phase structure for the dilatonic 
branes would be exactly the same as the non-dilatonic branes. The reduced Euclidean 
action 1% can be seen from (110) to take exactly the same form as the non-dilatonic branes 
given in (27), (28) and (29) in terms of the new parameters (118). 

So far, we have not addressed the issue regarding the validity of using the effective 
action in describing the phase structure of black p-branes throughout the parameter space 
considered 12 . Here we will address this for both non-dilatonic and dilatonic branes to- 
gether. For non-dilatonic branes, we need to keep the curvature of black brane spacetime 
uniformly weak throughout the parameter space. For dilatonic branes, in addition, we 



12 



We thank the anonymous referee for raising this concern. 
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also need to keep the effective string coupling uniformly weak. The effective string cou- 
pling can be read from (92) and the curvature can be calculated from the metric given in 
the same equation as 



R 



2 



a/2 



a^A_ 
4 A. 



2d - - 



P 2 



d — d 
D-2 



1? 



d+l 



(120) 



where g s = is the asymptotic string coupling. Note that for the scalar curvature each 
term in the square bracket is less than unity since p > r + > f_, A + /A_ < 1, a 2 /4 < 1 
and (d — d)/(D — 2) < 1. In other words, the square bracket contributes at most a factor 
of order unity to the curvature. So in order to keep the curvature uniformly weak, we 
need to have 

I 2 



l z R 



2/d 



«1, 



'121 



V2nQ d e a M 2 A a y\ 

where I is the relevant length scale under consideration, for example, it is the Planck scale 
l p in eleven dimensions or the string scale l s in ten dimensions. Note that the charge 
quantization gives y/2KQde a ^ ^ 2 /fij+i ~ Nl d with the integer N labeling the number of 
branes. So the uniformly weak curvature condition is 



NA 



a 2 / 4 



> 1. 



:i22i 



For M branes, a = and a weak curvature is the only requirement which can be satisfied 
when N 3> 1. For D3 branes, we need in addition a small g s . For those branes, it doesn't 
appear that there is any constraint on the parameter space we considered. For dilatonic 
branes, the condition (122) for weak curvature can easily be satisfied for non-extremal 
branes, i.e., r + > r_, for given large enough N since A_ is finite for p > r + . Now the 
effective string coupling as given in (120) can remain small for small g s when a > and 
can also remain small if g s is chosen to be small enough for a < 0. For this case, the 
parameter x used in the text falls in the range q < x < 1. If we consider extremal branes, 
i.e., r + = r_, we then have to limit to the range p > r + = r_ so that the curvature 
remains small. This can also give a small effective string coupling even when a < by 
the same token. Now we have the range of parameter space as q < x < 1. In other words, 
for dilatonic branes, we can consider the edge state only at one end, i.e., either at x = q 
or at x — 1, but not at both. 

The discussion presented in this paper corresponds to the single-scalar A = 4 black 
branes given in [12] (do not confuse this A with the discriminant notation used in the 
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main text of this paper). However, we have investigated the other A = 3, 2, 1 single scalar 
black branes discussed in [12] as well and found that the basic phase structure remains 
the same as that of A = 4 black branes even though the minimal value of d at which 
there exists a critical charge q c depends on the value of the respective A. 
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